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Abstract 

This paper presents new results for the modulus of families of walks on a graph—a discrete analog 
of the modulus of curve families due to Beurling and Ahlfors. Particular attention is paid to the 
dependence of the modulus on its parameters. Modulus is shown to generalize (and interpolate 
among) three important quantities in graph theory: shortest path, effective resistance, and max- 
flow or min-cut. 
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1. Introduction 

Beginning with an example, consider the simple, undirected graph shown in Figure Such a 
graph may be used to model, for example, a contact network for use in predicting the spread of 
disease, where each vertex indicates an individual in a population and each edge indicates that the 
two corresponding individuals have direct contact with each other. The disease is then assumed 
to spread throughout the network according to some stochastic process. An important question in 
such a model is the following. How easily can individual s transmit the disease to individual t? In 
a sense, we are asking for a distance-like measure between s and t. 

One obvious approach would be to define distance in the classical graph distance sense: the 
distance between s and t is the length of the shortest path connecting them. In this case, s and 
t would be at distance I from each other. For this particular application, however, this definition 
seems inappropriate. In graph distance, s is equidistant from u and t. However, any reasonable 
disease model would indicate that a disease originating in s is more likely to spread to t than to u; 
there are more avenues for transmission from s to t. 

Taking this shortcoming into account, we might try to define proximity through the number of 
distinct transmission pathways. That is, we could count the maximum number of pairwise edge- 
disjoint paths from s to t: 3 in this case. By Menger’s theorem, this is equivalent to the minimum 
number of edges that can be removed from the graph in order to separate s from t. This min-cut 
sense of proximity addresses the earlier concern with graph distance, since s is now three times 
closer to t than to u. However, this measure is also lacking; t is now equally close to both s and v, 
which again seems counterintuitive for the application. 
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A more recent description of distance that bridges the gap between these two extremes is the 
resistance distance, see for instance [12]. The resistance distance between s and t is defined as 
the effective resistance between the corresponding nodes in an electrical resistor network wherein 
the vertices of the graph correspond to junctions and the edges to resistors of unit resistance. In 
resistance distance, s and t are 0.454 units apart, s and u are 1 unit apart, and t and v are 1.918 
units apart. 

Although resistance distance does give a more satisfying answer for the contact network ap¬ 
plication, a number of interesting questions remain. For example, although the shortest path and 
min-cut distances have straightforward extensions to directed graphs, effective resistance does not. 
So, what is the natural analog of resistance distance for directed graphs? Or consider the following 
related type of question for the graph in Figure It is clear that s plays a crucial role in the 
transmission of disease between u and t. But, how might we quantify the role that t plays in 
transmission between s and v? In other words, how rich is the set of pathways connecting s to u 
through tl And how does this compare to the set of all pathways connecting s to u? Or for that 
matter, given any interesting collection of walks on the graph, how does one quantify the richness 
of the collection? 

A potential answer to these questions comes in the form of the modulus—a discrete analog of 
the modulus of curve families dating at least as far back as Beurling |3| and Ahlfors [T]. Although 
discrete versions of the modulus have been introduced before, the goal of the present work is to 
make explicit the connection to more traditional graph theoretic concepts: 

• We treat the modulus in a general framework, accommodating directed or undirected graphs, 
with or without edge weights (Section]^. Although this is not entirely new (see Section]^ 
for a description of related work), we do adopt a particular viewpoint that lends itself well 
to graph theoretic interpretations of the modulus. 

• After framing the modulus as a convex optimization problem, we derive and interpret a dual 
formulation of the modulus (Section]^. This is similar to the duality of the max-flow and 
min-cut problems, but much more general than has been done before. 

• Using the primal-dual formulation, we establish three connections to graph theory, namely 
that computing the modulus in certain circumstances is equivalent to computing either the 
graph distance, min-cut distance, or resistance distance (Section]^. Moreover, we prove that 
the modulus (which depends on a continuous parameter p) effectively interpolates among 
these three concepts (Section]^. 
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In Section the modulus is defined through real-valued densities on the edge set of the 
graph. We show that these densities generalize the concept of potential drop across resistors 
in the computation of resistance distance (see Theorem 4.2) and analyze the behavior of these 
densities (Section]^. 


2. The modulus 

2 .1. Graphs, walks, and families of walks 

In what follows, G is a graph with vertex set V and edge set E. The sets of vertices and 
edges are assumed finite, with n = \V\ vertices and m = \E\ edges. The graph may be directed 
or undirected. In the former case, E is represented as a subset of ordered pairs of vertices, in the 
latter as a subset of unordered pairs. We adopt the notation e = (x, y) to represent an edge in 
either case, with the understanding that (x, y) = (y, x) in the undirected case. For any function 
/ : Fi —)■ M and any edge e = (x,y) we will use the notation /(e) and f{x,y) interchangeably. 

The discussion in this paper is restricted to simple graphs, though there are natural extensions 
to multigraphs and to other more general objects. In general, G is assumed to be weighted, with 
weight function a : E ^ (0,oo). A graph is called unweighted if cj = 1. For the purposes of this 
paper, then, a graph is defined through the triple G = iV, E, a). 

Using the standard terminology, a walk on G is represented as 

• a hnite string of vertices V 1 V 2 V 3 ... Vr+i in V, with the property that {vi,Vi+i) G E for 
f = l, 2 ,...,ror, equivalently, 

• the corresponding string of edges 6162 ... e,., where e* = {vi, Uj+i) for i = 1 , 2 ,..., r. 

In either case, we require that r > 1, so that a walk is considered to traverse at least one edge in 
the graph. 

To each walk 7 = 6162 ... e,. is associated the graph length, £( 7 ) := r. That is, £( 7 ) counts the 
number of “hops” made in traversing 7 . More generally, given any edge density p : Fi —>■ M, the 
p-length of 7 , Ipi'j), is defined as 

r 

i=l 

The graph length of 7 coincides with the p-length for the particular choice p = 1. 

Throughout this paper, T represents an arbitrary nonempty family of walks on G. A particularly 
useful class of walk family is the connecting family, denoted F(s, t), of walks originating at a vertex 
s gV and terminating at a distinct vertex t gV \ {s}. 

2.2. Densities and admissibility 

Given a walk family F and an edge density p : Fi —)> M, the graph length and p-length of F are 
defined as 

£{T) := inf £( 7 ) and £p{r) := (1) 

7Gr 7Gr 

respectively. 
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The family T determines three sets of densities called the admissible set ^(T), the relaxed 
admissible set A'(r), and the restricted admissible set ^*(r) dehned as follows. 


A'(r) 

A(T) 

A*(r) 


{p : .EM : £p(r) > 1, 0 < yo}, 

{yo : M : £p(r) > 1, 0<yO<l}. 


Note that these dehnitions are made in the order of decreasing size: ^*(r) C ^(T) C ^^(T). The 
reason for defining three different admissibility sets will be apparent by Remark |2.1[ 


2.3. Energy and modulus 

Given a real parameter p > 1 or p = oo, the p-energy of a density p is 

f if 1 < p < oo 

max |p(e)| if p = oo 

e&E 

The definition for is consistent in the sense that 

Vp ; S —>• M lim Ep{pY^^ = Eoo{p). 

p^oo 

For 1 < p < oo, the p-modulus of T is defined as 

Modp(r) := inf E„{p). 

^ peA(r) ^ 

(The modulus of an empty family Modp(0) is defined to be zero, since the choice p = 0 is trivially 
admissible.) 


2 . 4 . Extremal densities 

An edge density p G A(r) is called extremal for a given family T and a given p if 

Modp(r) = Ep{p). 

The notation p* or pp is frequently used to denote an extremal density, with the latter used in 
order to make the dependence on p explicit. 

Lemma 2.1. Given T and 1 < p < 00 , there exists an extremal density p*. //1 < p < 00 , then p* 
is unique. 

Proof. For any 1 < p < 00 , the function p 1 —)• Ep{pY/^ is a norm on the m-dimensional space of 
edge densities, as is the function p 1 —)■ foo(p)- Moreover, the set A(r) is closed and convex. Thus, 
the modulus problem can be restated as the problem of finding a point in a closed convex subset 
of M'” that is closest to the origin in a particular norm. Such a point always exists, and is unique 
provided the norm is strictly convex, as it is in the cases 1 < p < 00 . □ 

The following observation, which shows that the nonnegativity constraint on p need not be 
explicitly enforced, will prove useful in Section 
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Lemma 2.2. Let F and 1 < p < oo be given. Then 


Modp(r) = inf Snip)- 
^ peA'(r) ^ 


( 2 ) 


An extremal p* G ^^(F) always exists for the relaxed problem, and is unique if 1 < p < oo. 
Moreover, for 1 < p < oo any such extremal p* satisfies the bounds 0 < p* < 1, and thus p* G A*(F). 
When p = oo there exists at least one extremal p* G 74 *(F). 

Proof. Existence and uniqueness (for 1 < p < oo) of a minimizer for the infimum in Equation Q 
are proved as in Lemma |2.1[ 

Now, suppose p G ^^(F) takes a negative value on some edge e' G E, and define p'^ as p'^{e) = 
max{/ 9 (e), 0}. Since ^p+{'y) > £p(j) > 1 for any 7 G F, /?+ is also admissible. Moreover, on every 
edge e € E, |/ 9 ^(e)| < \p{e)\ and \p'^{e')\ = 0 < \p{e')\. For 1 < p < 00 , this implies £p{p'^) < £p{p), 
showing that any minimizing p* is in ^(F). For the p = 00 case, the inequality Sooip~^) ^ £oo{p) 
holds, but need not be strict. However, a minimizing p* can always be found in H(F). 

Next, suppose p G H(F) and that p{e') > 1 on some edge e' G E. Define pfe) = m.m.{p{e), 1}. 
For any 7 G F, either = f’p( 7 ) > 1, or 7 traverses an edge where takes the value 1. In the 

latter case, the inequality ipf'j) > 1 is trivial, so p' G H(F). Since p'{e) < p{e) on all edges, and 
since the inequality is strict on at least one edge, we have £p{p') < £p{p) (for 1 < p < 00 ). This 
proves that p* <1. □ 


Remark 2.1. A consequence of Lemma 2.2 is that the admissibility set A(F) in the definition of 
modulus can be replaced with either a more relaxed set, A'(F), or a more restricted set, A*(F), 
without changing the outcome. That is. 


Modp(F) := inf Spip) = inf Spip) = inf Spip) 
^ peA(r) ^ peA'(r) ^ peA*(r) ^ 


for all 1 < p < 00 . 


A convenient choice of admissibility set can simplify proofs in certain circumstances. For example, 
when deriving the Lagrange dual problem in the case 1 < p < 00 (Section it is better to use 
the most relaxed set A'(F). The proof of Theorem 3.1 utilizes the lower bound in A(F), while the 


bounds of Section make use of both upper and lower bounds in A* (F). 

2.5. An instructive example 

As described in the introduction, the p-modulus is a generalization of three fundamental quan¬ 
tities in graph theory: shortest path length, effective conductance, and min-cut. Before proving 
these facts, it is useful to explore an example. Consider the (unweighted and undirected) graph 
shown in Figure [2]— comprising k parallel simple paths of i hops connecting node s to node t —with 


F = F(s,t), the connecting family of walks defined at the end of Section 2.1 

It is a straightforward exercise to show that the choice p* = j extremal for all 1 < p < 00 . 
Thus, for 1 < p < 00 , the modulus of F is 


Modp(F)=Tp(p*) = A:£(^y=^, 


and for p = 00 , 


Modoo(F) =Too(p*) = j. 
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Figure 2: A graph consisting of k simple paths with I hops each, connecting s to t. 


There are several interesting observations to be made. First, as a function of p, the modulus is 
continuous, monotone decreasing, and decays to zero like as p —)• oo if .£ > 1. All three of these 
properties are generic, as shown in Section 

Moreover, for 1 < p < oo, the modulus depends on both the number of distinct walks, k, and 
on the length of these walks, i. The modulus is large for F containing many, short walks, and small 
for F containing few, long walks. For p near 1, the modulus is much more sensitive to changes in 
k than to changes in while, for p very large, the modulus is much more sensitive to changes in 
i. In the extreme case p = 1, the modulus loses all dependence on l\ a large family of walks has 
large modulus regardless of their lengths. On the other hand, when p = oo the modulus loses all 
dependence on k: a family’s modulus depends only on the length of the shortest walk and not at 
all on the number of walks. The case p = 2 is also interesting. In this case, the modulus is k/i, 
which is the effective conductance of the graph viewed as a resistor network between s and t with 
unit resistors placed on each edge. 

Again, these are general properties. Section proves that, in a specific sense, the 1-modulus is 
a generalization of the min-cut problem, that the oo-modulus is always equal to the reciprocal of 
the length of the shortest walk, and that the 2-modulus is a generalization of effective conductance. 
The results of Sectionthen, show that p can be thought of as a tuning parameter, which adjusts 
the balance of importance of “many walks” with “short walks” in the modulus. 

3. Lagrange dual formulation 

The p-modulus problem is a convex optimization problem: it involves minimizing the convex 
function £p over the convex set A(r) defined by a set of linear inequalities. Furthermore, it was 
shown in [3] that this set of inequalities can be assumed finite; more precisely, the following theorem 
is true. 

Theorem 3.1 ([3])- Let T be a given family of walks on a graph. There exists a finite subfamily 
F* C F such that A{r*) = A(r). 

Such a finite subfamily is called an essential subfamily for the modulus problem. Ensuring that 
^pil) > 1 is true for all walks in the essential subfamily is equivalent to ensuring the inequality 
for the entire family. Thus, even when a described set of walks (e.g., the connecting family T{s,t)) 
contains infinitely many members, it can always be assumed that F has been replaced by a finite 
essential subfamily (e.g., the simple paths from s to t when F = r(s,t)). With this assumption. 
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the modulus problem becomes an ordinary convex program, in the sense of m Sec. 28], taking 
the form 


minimize £p{p) 

subject to J\f , e)p{e) >1 V 7 G F, 


(3) 


eeE 

p{e) >0 Ve G FI, 

where AA( 7 , e) is the |r| x \E\ matrix whose ( 7 ,e) entry contains the number of times walk 7 
traverses edge e. 

Using standard techniques from convex optimization, we define the Lagrangian for the problem 
as 


Lp{p,X,p) := £p{p) + M - ^A/'(7,e)/9(e) j -^p{e)p{e), 

7er V eSE / eSE 


(4) 


where A : F —)• [0, 00 ) and p : E ^ [0, 00 ) are the Lagrange dual variables associated with the 
admissibility constraints. The usual weak duality statement takes the following form. 

Lemma 3.2. Let F and p be given. Then 

Modp(F) = inf sup Lp{p,X)> sup inf Lp{p,X)= sup Ep{X, p) 

P-E^R X-.r^[0,cxi) A:r->-[0,oo) A:r^.[0,oo) 

/f:E^[0,oo) >-[0,oo) /f:E^[0,oo) 


where 


Fp{X,p)-.= ini Lp{p,X,p). 

p-.E^R 


Given any A : F —>■ [0,oo) and any p : E ^ [0,oo), weak duality implies that Modp(F) > 
^^^(A, p). The problem of maximizing this lower bound over all nonnegative A and p is the Lagrange 
dual problem to the modulus: 

maximize Ep{X, p) 
subject to A( 7 ) >0 V 7 G F, 
p{e) >0 Ve G FI. 

Remark 3.1. Because the inequality constraints in ^ are affine, it follows from Theorem 28.2 
of [H] that strong duality holds. That is, there exists a saddle point {p*, X*, p*) of the Lagrangian. 
This is an important property of the modulus problem, especially in the context of numerical 
approximation, and was explored in a specialized setting in [3]. For the present paper, the weak 
duality result is usually sufficient. 

When working with the dual problem, it is convenient to have a more explicit formulation. We 
derive such formulations in three different cases: l<p<oo,p = l and p = 00 . 

3.1. The case 1 < p < 00 

When 1 < p < 00 , Lemma [2 . 2| implies that the p > 0 constraints need not be explicitly enforced 
in Q, and can thus be removed. The simplified Lagrangian (omitting the dual variables p as they 
are no longer necessary) takes the form 

Lp{p^^) ■= + ^ “ '^^{XF)p{e)\ , 

eSE 7er V eSE / 

7 


which is differentiable in each p{e), since p > 1. Thus, the value of the dual energy 


Jp(A) = inf^l ^cj(e)|p(e)|P + ^A(7) - ^7^(7,e)p(e) j > 

^ 7er V eS-B / 

can be found by solving VpLp = 0. The optimal density, p\, for a given A must satisfy the 
stationarity condition 


1 

p-i 


pa{e)px{e)\px{e)Y' ^ ^AA(7,e)A(7) = 0 \/eeE. 

7er 

Since A > 0 and J\f has non-negative entries, the above equation implies that pa > 0 and so 

Pxie) = 

After some simplification, this produces the following formula for Ep{\). 

Jp(A) = L{px, A) = ^ A( 7 ) - (p - 1 ) ^ a(e) j e)A( 7 ) 


p 

p-i 




7 Gr 


We conclude that the Lagrange dual optimization problem can be written as 

1 


p 

p-i 


maximize ^A( 7 ) - (p- 1 ) ^cr(e) 
7er eSB 

subject to A( 7 ) >0 V 7 G T. 


pa{e] 


^AA(7,e)A(7) 


7 er 


(5) 


( 6 ) 


Remark 3.2. Because strong duality holds (see Remark 3.1), solving ^ is equivalent to computing 
the modulus. This provides an interesting interpretation of the modulus: rather than assigning 
values to edges (through p), we instead assign values to walks in T (through A) in such a way as to 
maximize the dual energy. Equation ([^ provides the formula for determining the extremal density 
p* from an extremal A*. On the other hand, the extremal A* are typically not unique. 


3.2. The case p = 1 

In the case p = 1, the general Lagrangian in Equation Q is not differentiable, so the preceding 
argument does not apply. This issue can be circumvented by replacing the 1-modulus optimization 
problem with an equivalent problem: 

minimize E o-(e)p(e) 

eSB 

subject to AA( 7 , e)p(e) >1 V 7 G T 

eSB 

p(e) >0 Ve G E. 
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Dropping the absolute values from the energy has no effect on the outcome, since the values of p are 
required to be nonnegative. However, without the absolute values, the energy is now differentiable 
for general densities p : i? —)■ M. 

Introducing the dual variables A : F —)■ [0,oo) and p : E ^ [0,oo), the Lagrangian for this 
problem is 


Li{p,X,p) := '^a{e)p{e) + ^A( 7 ) | 1 - J\f {-f, e) p{e)\ - '^p{e)p{e) 

e&E 7er V e&E ) e^E 

= + fo-(e) - ^AA(7,e)A(7) - ^(e)^ . 


7 er 


e&E 


7 er 


Since the Lagrangian is affine in p(e), it is clear that 


E ^(7) if crie) - Af{7,e)X{-f) - pie) = 0 \fe e E 

Ei{X,p) := inf ^Liip,X, p) = { 76 r 7 er 

—00 otherwise 


p:E^R 


The dual problem involves maximizing Ei over feasible A and p, and therefore it will never be 
advantageous to choose the dual variables such that T'i(A, p) = —00 (assuming this can be avoided). 
In fact, provided A satisfies the inequalities E 7 er-^( 7 T)^( 7 ) < <7(e) (e.g., A = 0), a suitable 
nonnegative p can always be defined to make T'i(A, p) hnite, whereas, if any of the aforementioned 
inequalities is violated in the choice of A, there will be no nonnegative p for which Ei{X, p) does 
not evaluate to — 00 . 

Thus, the dual problem can be written as follows. 

maximize 

7 er 

subject to E AA( 7 , e)A( 7 ) < (T(e) \/e e E (7) 

7 er 

A( 7 ) >0 V 7 G T. 


3.3. The case p = 00 

In the p = 00 case, the general Lagrangian is once again not differentiable. In this case, an 
appropriate equivalent problem is the following. 

minimize t 

subject to Mi'j, e)pie) >1 V 7 G T 

e£E 

0 < pie) <t \/e € E. 

The associated Lagrangian takes the form 


Loo(L p, A, p,-!]) :=t + J2 ^( 7 ) 
7 er 


1 - ^A7(7,e)/)(e) 

e£E 


eeE e£E 



- y^-^(7, e)A(7) - pie) + vie) 


7 er 


7 er 
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As in the p = 1 case, the Lagrangian is affine in the primal variables t and p, so the dual 
objective function Joo(A,/r,r/) is given by 

f E A( 7 ) if E ^(e) = 1 and A^( 7 , e)A( 7 ) +/i(e) - 77 (e) = 0 Ve G ^ 

J'oo(A,/U,r7) := < 76r esE 7er 

1—00 otherwise 

Similarly to the case p = 1, this produces the following form of the dual problem. 

maximize 5ZA(7) 

7 er 

subject to ^AA( 7 , e)A( 7 ) < 77(e) \/e € E 

7 er ( 8 ) 

^ 77(e) = 1 

e£E 

a( 7) >0 V7 G r. 


4. Connections to common graph theoretic qnantities 

One of the key contributions of the present work is to make explicit the connections between 
the concept of p-modulus and the graph theoretic concepts of shortest path, effective conductance, 
and min-cut. The present section is devoted to establishing these connections. 


4-1- Shortest path 

The connection to shortest path is the broadest and simplest result to establish. Recall that 
£(r), defined in Equation is exactly the length of the shortest walk in T (in terms of hops). 

Theorem 4.1. LetT he any family of walks, then 


Modoo(r) 


1 


Proof. Let 7 ' G T be a walk such that ^( 7 ') = ^(T) = £. The density po = 1/i is admissible, since 
for any 7 G T 

^po(7) = '^J^{l^e)pQ{e) = > 1. 

eeE ' 

Thus Modoo(r) < £~^. On the other hand, if p G A(r), then 


1 < £p(7') = Y^-^( 7 ',e)p(e) < £( 7 ') maxp(e), 

e£E 

SO every admissible density must satisfy £oo{p) > £~^- 


(9) 

□ 
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4-2. Effective conductance 

An undirected graph G = {V,E,a) is a model for a resistor network with edges representing 
resistors with conductances a, connected at junctions represented by the vertices. Let s and t be 
distinct vertices in V. The effective conductance Ceg(s, t) (the reciprocal of the effective resistance) 
between s and t can be found by minimizing the total power 

Power = y) ('^(^) - 

{x,y)£E 

over all voltage potentials <() : P —>■ M satisfying 4>{s) = 0 and (f){t) = 1. The following theorem 
shows that, for 1 < p < oo, the extremal density p* for the modulus of a connecting family of walks 
can be related to a generalized voltage potential. Such a result, in the language of extremal distance 
and in the special case p = 2, first appeared in the work of Duffin [6] (see also [71 Proposition 6.2]). 
For a version in metric spaces, see m Theorem 7.31]. 

Theorem 4.2. Let G = {V,E,a) be an undirected graph, let P = r(s,f) he the conneeting family 
of walks from s to t, two distinct vertices in V, and let 1 < p < oo. Let p* be the extremal density 
for Equation (©• Then there exists a vertex potential (/>* : P —)■ M such that (j)*{s) = 0, </>*(!) = 1, 
and 

p*{x,y) = \(l)*{x) - cj)*{y)\ y{x,y)£E. (10) 

Moreover^ this solves the optimization problem 

minimize E a{x,y)\4>{x) - (j){y)\P 

ix,y)&E 

subject to (j){s) = 0, 4 >{t) = 1. 


When p = 2, it follows that Mod 2 (r(s,f)) = (765(5, f). 

Proof. Let p* G ^(T) be the extremal density. Note that lp*{T) = 1, for, if not, the density 
p' = p* ftp* iV) is also admissible and has lower p-energy than p*. Define (f*{s) = 0 and for 
X G P \ {s} 

4 >*{x) = min 
7er(s,x) 


Since P = r{s,t), (j)*{t) = 1. To see that Equation (10) holds, define p'{x,y) = \(j)*{x) — 4>*iy)\ for 
(x, y) G E. For any 7 = SV 2 V 3 ... Vrt G P, 


V(7) = l</>*('S) - f*iv2)\ + \f}*iv2) - (Tivsf H-h \(l)*{Vr) - f)* (tf > \f*{t) - ())*(5)| = 1, 


SO p' G A(r). Let (x,y) G E be an edge. Without loss of generality, <f)*{x) < 4>*{y) and y s. We 
claim that 

0 < p'{x,y) = (Tiy) - (/>*(x) < p*{x,y). (11) 

The case x = s is trivial. Suppose that x s and let 7 G P(s, x) be a walk such that ip*{'y) = 4>*{x). 
Letting 7 ' G P(5, y) be the walk obtained by first traversing 7 and then traversing the edge (x, y), 
we have 

4 >* {y) < V (7O = V (7) + p*{x,y), 


which implies Equation (11). 
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Equation ( [IT| ) implies that £p{p') < £p{p*) and, by uniqueness of the extremal density, p* = p'. 
To see that (jf solves the optimization problem, let 0' : E —)• M with (t)'{s) = 0 and ^'{t) = 1 and 
define p'{x,y) = \(p'{x) — As before, p' G A(r) and so 

^ ^pip') = ^ (^{x,y)\(t>'{x) - 4>\y)\p. 

{x,y)&E {x,y)&E 


□ 


4-3. Max-Flow Min-Cut 

Modulus is also closely related to the Max-Flow Min-Cut Theorem [ 8 ]. There are a number of 
different ways to see this; we shall focus on the connection to Ford and Fulkerson’s original work. 
Let G = (E, E, a) be a weighted, undirected graph and let s, t G E be distinct vertices. Let T 
be the family of simple paths from s to t. (T C r(s,t) is an essential subfamily in the sense of 
Theorem 3.1 ) In the terminology of Ford and Fulkerson, the positive edge function a is called the 
capacity of the edge. A flow can be thought of as a function A : T —)■ [0, oo) with the property that 


^ AA( 7 , e)A( 7 ) < cr(e) Ve G E, 
7 er 


and the value of such a flow the sum of A( 7 ) over all 7 G T. Thus, Ford and Fulkerson’s maximal 
flow problem is equivalent to Equation Q and forms a lower bound for the 1-modulus of T. 

A disconnecting set is a subset of edges D <Z E with the property that each 7 G T traverses 
at least one edge in D, and a cut is a disconnecting set that contains no other disconnecting sets 
as proper subsets. The value of a disconnecting set D, denoted v{D), is defined as the sum of the 
capacities of all edges in D. The Max-Flow Min-Cut Theorem is stated as follows. 

Theorem 4.3 (Max-Flow Min-Cut). The values of the maximal flow and minimal eut on any 
graph are equal. 


The connection between this and the 1-modulus is established in the following theorem. 

Theorem 4.4. Let G = (E, E, a) be a weighted, undirected graph, let s,t G V be two distinct 
vertices and let T{s,t) be the conneeting family of walks between s and t. Then Modi(r(s,t)) is 
equal to the value of maximum flow (and minimum cut) from s to t. 

Proof. Let T C r(s,t) be the family of simple paths connecting s to t. T is an essential subfamily 
and so Modi(r(s,t)) = Modi(r), which is bounded below by the value of the maximum flow, as 
we saw above. 

Let D E E he a disconnecting set of minimal value v{D) and define the edge density p as 


P(e) 


1 if e G D 
0 otherwise 


Since D \s a disconnecting set, p G A(r(s,t)). Moreover, 

Modi(r(s,t)) < £i{p) = Y^cj{e)p{e) = Y,cj{e) = v{D), 

e£E e£D 

providing the opposite inequality. □ 
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5. Dependence on the parameter p 


The results of the preceding section are that the p-modulus can be thought of as a generalization 
of shortest path {p = oo), effective conductance {p = 2), and max-flow or min-cut {p = 1). Another 
way to see this is to consider an undirected, unweighted graph G and let T = r(s,t) be the 


connecting family of walks between distinct vertices s and t. (Recall the example in Section 2.5 ) 


In this setting, Modoo(r) equals to the reciprocal of the length of the shortest path in T, while 
Modi(r) equals to the maximum number of edge-disjoint paths from s to t (a consequence of 


Theorem 4.4 and Menger’s theorem). Mod 2 (r) strikes a balance between the other two, in a sense 


preferring both shortness and abundance in the collection of paths connecting s to t. 

In fact, as the present section demonstrates, the parameter p can be thought of as a continuous 
“tuning parameter”, which balances the sensitivity of the modulus between walk length and walk 
abundance. More precisely, given any graph G = {V,E,a) and any nonempty family of walks T, the 
functionp i—>■ Modp(r) is continuous and decreasing. Moreover, Modp(r)^/P converges to Modoo(r) 
as p ^ oo. In fact, with the correct normalization, the convergence is monotone (see Theorem |5.2[ ). 
We begin with bounds on the modulus. 

Lemma 5.1. Let G = {V, E, a) he a graph with weights a > 0. Let the family T and the parameter 
1 < p < oo be given. Then 

^ < Mod,(r) < 

where £ := i{T), 

o’min ■= min(T(e), and criE) := > c(e). 
eeE 

e£E 

Proof. The upper bound can be obtained by considering the edge density p = This density is 
admissible, because for all 7 G T, ipi'j) = ^£( 7 ) > 1- Thus 


Modp(r) < Lpip) = Y^(t( 


eeE 


^ £P 


4^ 

£p 


The lower bound can be obtained as follows. Let 7 ^ G T such that £(j^) = £(r) = £. By 
Equation ([^, for every p G ^(T) there exists an edge e G E such that p{e) > £~^. Thus, 


^P^p) = / )cr(e)p(e)^ > O-minp(e)^ > 


o-n 


e&E 


£P 


□ 


Theorem 5.2. Let G = {V,E,a) andV he given. Then 

lim Modp(r)p = = Modoo(r). 

p^oo £\r) 

Moreover, the following monotonicity properties holds for 1 < p < q < 00 : 

Modg(r) < Modp(r), 

a(E)-i/PModp(r)Vp < fT(E)-i/'?Modg(r)i/'?. 
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( 12 ) 


(13) 

(14) 









Remark 5.1. Observe that (12) and ( |14[ ) combine to yield 

a{E)-^/P ModpiVf/P t Modoo(r) as p —)• oo. 


( 15 ) 


Proof. Taking the pth root of the inequalities in Lemma 5.1 proves (12). Equation (13) follows 
from the fact that for any p G 74*(r) we have 0 < p < 1, and since q > p, this implies that 


£qip) = ^{e)p{ey < Y ^(e)/3(e)^ = £pip)- 

e£E e£E 


(16) 


To prove monotonicity in (14) we use Holder’s inequality. Given p G A*(r), assume that 1 < p < 
q < oo. Then, using the conjugate exponents q/p and q/{q — p), 


£p{p) = Y = Y ^(e)"/V(e)" • a(e)i-^/'? 

eeE eeE 

/ \plq / \ i-p/g 

< ^ <T(e)p(e)’ ^ a(e) = 


(17) 


\e^E 


\eeE 


This can be rewritten as 


a{E)-^/PSp{py/P < u(E)-V'?T,(p)Vt 
Picking p to be extremal for Modq(r) we get 

Theorem 5.3. Let G = (E, E, a) and T be given. Then the the function 

p I—)• Modp(r) 


(18) 

□ 


is continuous for 1 < p < oo. 


Proof. If we fix p and let q ip, from (13) we get 


limsupModq(r) < Modp(r), 

qip 


and from (14) we get 


liminf Modg(r) > liminf ct(E)^ Modp(r)'^/^ = Modp(r). 

qlp qlp 


So limgj^pModq(r) = Modp(r). 

Likewise, ii qf p, then from (13) we get 


liminf Modq(r) > Modp(r), 


and from (14) we get 


limsupModg(r) < limsupcr(E)^ Modp(r)'^/^ = Modp(r). 

<7tp qtp 


So limg_^p Modq(r) = Modp(r). 


□ 


14 











6. Properties of extremal densities 


We end with two properties of the extremal density in modulus problems: it is continuous as a 
function of p, and it is closely related to the gradient of the modulus with respect to the weights 

a. 


6.1. Continuity of the extremal density 

This continuity result arises from Clarkson’s inequalities. (See Theorem 2 of [5].) For a given 
p > 1, let II • lip represent the standard p-norm on M™, defined on a density p : E as 

II* (E I*)!" 

\eeE 



In this context, the relevant Clarkson inequalities take the following form. 

Theorem 6.1. Let p > 1, p' = pj{j> — 1), x, y G M™, and || • ||p be the standard p-norm. 
If P>‘^, then 


x + y 
2 


p 

+ 

p 




p 

p- 


If 1 < p < 2, then 


x + y 


+ 


X -y 


< (^ll^llp + ^l 


A first consequence of Clarkson’s inequalities is that admissible densities that are almost- 
minimizers of the energy are close to the extremal density in p-norm. 


Lemma 6.2. Let G = {V,E,a) he a graph with 0 < U min < a, and let T be a family of walks on 
G. Assume p > 1, and let p* he the unique extremal density for the p-modulus ofT. Then, given 
any admissible density p G A(r) the following holds. 

If P>‘I, then 

2P-1 

Up - p* % < — {£p{p) - Modp(r)). (19) 

< 7 min 

If 1 < p <2, then 


\\P-P*\% < 


2 P 


(Jr, 


fModp(r) + ^(£:p(p)-Modp(r)) 


— Modp(r) p 


( 20 ) 


Proof. It is convenient to define two additional edge densities, / and /*, as 

f{e):=a{e)pp{e) and f*{e) := a{e)pp*(e). 

Then 


P = Sp{p) and \\r\\P = Sp{p*) = ModpiT), 


and 


f-f* 


' = i E - p’leT >^\\p- P'K 

P e&E 
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( 21 ) 

( 22 ) 
































Also, 


f + r ^ 

2 P 


eS-E 


p{e)+p*{e) 


— f 

- Or) 


/ 9 (e) + p*(e) 


> Modp(r), 


(23) 


where the final inequality follows from the fact that the admissible set A(r) is convex, so that 

(p + p*)/2 G ^(r). 

When p > 2, Equations (21) through (23) combined with Theorem 6.1 show that 

<^min I 


which proves 


2 P ^\\p - p* II? < + 2 Modp(r) - Modp(r), 

When 1 < p <2, a similar argument shows that 


a 


p_ 

min \ p 


2 P 


\\p-p*rp < -£:p(/o) + -Modp(r) -Modp(r)T 


implying (20). 


□ 


Theorem 6.3. Let T be a family of walks on a graph G and let 1 < p < oo. Then 

lim||pp-pq|| =0, 
q^p 

where pp and pg are the unique extremal densities for p and q > 1 respectively. 

Note that the norm in the theorem is irrelevant, since all norms on are equivalent. 


Proof. By Lemma 6.2, it is enough to show that £p{pq) —)• Modp(r) as g —)• p. For q > p, 
Equation ( [l7| ) shows that 

Modp(r) < Spipg) < a{E)^--^£g{pg)-^ = a(E)^-f Modg(r)? 


and for q < p, Equation 


'P\HqJ 
shows that 


Modp(r) < £p{pq) < £q{Pq) = Modg(r). 


Then, by continuity (Theorem 5.3), letting q tend to p yields the result. 


□ 


6.2. Extremal density as a gradient 

We now turn our attention to the behavior of the modulus as a function of the weights a : 
E —7> (0,oo). For this discussion, we fix a graph G = {V,E) and a non-empty family of walks T, 
but allow the edge weights, a, to vary. In order to make explicit the dependence on the weights a, 
we use the notation £p{p] a) and Modp(r; a) to denote the p-energy of p and the p-modulus of T 
respectively. We restrict attention to p < oo, since £^0 and Modoo do not depend on a. The next 
two lemmas show that modulus is continuous and concave in a. 


Lemma 6.4. Let 1 < p < 00 . The function a 1 —)• Modp(r;(T) is Lipschitz continuous. 
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Proof. Let cti,cj 2 : E —>• (0,cx)) and let pi and p 2 be extremal densities for Modp(r;(Ti) and 
Modp(r; 1 T 2 ) respectively. Without loss of generality, assume Modp(r;cJi) < Modp(r;(T 2 ). From 
Lemma 2.2, it follows that \pi\ < 1, so 

Modp(r;cT2) - Modp(r;cTi) < £p{pi](T2) - £p{pi]ai) = ^ cT2(e)pi(e)^ - '^ai{e)pi{eY 

e£E e£E 

= ^(0-2(6) - ai{e))pi{ey < ^ |fT2(e) - cJi(e)| < m||cj2 - cti||oo, 

eSE eSE 


where 


is the standard max norm on 


□ 


Lemma 6 . 5 . Let 1 <p < 00. The function a 1—)• Modp(r;(T) is concave. 

Proof. Let ai,a2 ■ E ^ (0,oo) and let 9 G [0,1]. Let pi, p2 and p* be extremal densities for cri, 02 
and a* := 9ai + (1 — 9 )a2, respectively. Then 

0Modp(r;cJi) + (1 -0)Modp(r;fT2) < 9 {e)p*{eY + (1 - 0) ^ cf2{e)p*{eY 

eSE eSE 

= ^ cr*(e)/3*(e)^ = Modp(r;6*(Ti + (1 - 6*)cr2). 

eSE 

□ 


For p > 1, the following theorem provides an interpretation of the unique extremal density in 
terms of the gradient of Modp(r; a) with respect to cr. 

Theorem 6.6. Let 1 < p < 00 and define the function F(a) := Modp(r;(T) on weights a : E ^ 
( 0 , 00 ). Given a weight a, let p* denote the unigue extremal density for Modp(r;(T). Then 


dEja) 

da{e) 


Plief 


(24) 


Proof. Let cr : Fi —>• (0, 00 ) and r] : E ^ M. For h > 0 sufficiently small, ah ■= a + hp is positive. Let 
be the extremal density for Modp(r; ah), and let p* = Pq, the extremal density for Modp(r; a). 
Note that, for any p : Fi —)• M, 


Spip', (^h) = + hp{e))\p{e)\P = £p{p; u) + F ^ p{e)\p{e)\P. 

eSE eSE 


In particular 


Modp(r;cj,j) < £p{p*-,ah) = £p{p*-,cr) + h'^p{e)p*{e)P = Modp(r;c7) + F ^ p(e)p*(e)^’ 

eSE eSE 


and 


Modp(r; cT/j) = £p{p*h,ah) = £p{pl-,a) + h'^p{e)p*h{e)P > Modp(r; u) + F ^ p(e)p^(e)^’. 

eSE eSE 

Combined, these two inequalities imply that for F > 0 sufficiently small 

E s <Y^,(e)p- (er. (25) 

eSE eSE 
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Now, since each p*^ G A*(r) by Lemma 2.2 


£p{ph - Modp(r ; fj) = Sp{pl; ^ p{e)pl{ef - Modp(r; a) 

eGE 

< I Modp{T;ah) - Modp(r;cj)| +h'^\p{e)\. 

eGE 


Lemma 6.4 implies that the right-hand side vanishes in the limit as h —)• 0, so £p{p’^',a) 
Modp(r;iT). And Lemma 6.2 then implies that p\^ p*. Taking the limit in (25) yields 




h^0+ 


h 


eGE 


which, since p was arbitrary, implies (24) 


□ 


7. Related work 

The concept of the modulus of a family of walks on a graph is a natural analog of the continuum 
dehnition of modulus of a family of curves in the plane, introduced by Beurling [1] and Ahlfors [I]. 
The task of detailing its historical development is complicated by the facts that modulus seems 
almost infinitely generalizable and that similar ideas exist in different fields, with different notation 
and terminology. Certainly, the work on maximal flow by Ford and Fulkerson [8] is closely related 


to modulus (see Theorem 4.4), as is the work on extremal length in networks by Duffin [6]. Indeed, 
Dufhn’s paper cites Beurling and Ahlfors as well as Ford and Fulkerson in the introduction, and 
Dufhn’s extremal width is equivalent to a particular modulus problem. (Compare Theorem 2 of [B] 


to Theorem 4.2 of the present work.) More recently, one finds the work of Schramm which 
uses a form of modulus to study square tilings, and the work of Haissinsky cni, which develops 
a discrete version of the modulus quite similar to the dehnition in the present work, but with a 
different goal in mind. 


8. Discussion 

This paper presents an exploration, in a fairly general setting, of the ways in which modulus 
generalizes and expands upon important graph theoretical concepts. Modulus is shown to generalize 
three fundamental concepts from graph theory—shortest path, effective conductance, and min- 
cut—and to interpolate, in a sense, among these three concepts. These results depend upon lower 
and upper bounds on the modulus, which are established using tools from convex optimization. 

With these connections hrmly in place, a number of questions naturally arise. For example, 
resistance distance has proven quite valuable in a wide variety of applications. As we have shown, 
on a weighted, undirected graph, this distance is equivalent to the 2-modulus, in the sense that the 
resistance distance between s and t in 1/ is Mod 2 (F(x, y))“^. Is p = 2 the only useful case? Or are 
there applications for which some other value of p is useful? Moreover, there is no obvious sense of 
effective resistance on directed graphs, but such a concept may well be useful (see, e.g., [IS]-) We 
are aware, for example, of a recent arXiv paper m which constructs a notion of resistance distance 
on directed graphs that preserves certain important control theoretic properties. The average of 
Modp(F(x, y)) and Modp(F(y, x)) might also be an interesting quantity to look at in applications. 
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There are also several interesting open questions regarding the dependence of the modulus and 
extremal density on p beyond the continuity and monotonicity proved in Section For example, 
it is not yet clear how smooth the function p i—)• Modp(r) is in general. Moreover, it would be 
interesting to better understand the behavior of the unique extremal density pp either as p —>■ 1 or 
as p —7- oo. In both cases, one might wonder if there is a limiting density and, if so, if the limiting 
density is extremal for the 1- or oo-modulus. 

Another potentially interesting direction is graph optimization, similar in spirit to the results 
of Reference [9]. As shown in Section 6.2, the function a i—)• Modp(r; a) is concave, so the problem 


of maximizing the modulus subject to convex constraints on a (e.g., that the sum of a over all 
edges is bounded) is a convex optimization problem whose solution could give some insight into 
the structure of a graph relative to a given family of walks. 

As mentioned in the previous section, generalizations of the modulus presented here seem almost 
endless. Some obvious directions of research include the replacement of edge densities with vertex 
densities, or using a mixture of both. Moreover, the objects in F don’t actually need to be walks. 
One might consider a “cut modulus”, for example, where F is the family of all edge cuts separating 
two given sets of vertices. 

Another interesting line of inquiry lies in understanding what the modulus can say about the 
structure of a graph. As mentioned in the introduction, for example, it can be interesting to 
consider moduli of families of walks from s to t passing through u. One might also consider other 
families, such as the family of simple paths with at least L hops, or the family of simple cycles 
with diameter at least D. What might the moduli of these families say about the structure of the 
graph? What is the expected modulus of these families on a given random class of graph? 

Perhaps even more intriguing is the interpretation offered by the dual formulation of Section 
In the dual interpretation, the edge density p is replaced by non-negative function A on F. Through 
duality, A is related to the sensitivity of the modulus to the admissibility constraints, giving a 
natural ranking of the importance of the walks in F. Those walks with larger A are more important 
in the modulus problem than those with smaller A. This perspective of most important walks is 
rich for exploration and has been developed in [2]. For instance, by replacing F by some properly 
chosen subset F' with large A values, one can obtain a subfamily with nearly the same modulus as 
F. This viewpoint also lends itself to a concept of sparsification subordinate to a modulus problem: 
by choosing a subfamily F' of most important walks, and then removing from E all edges that do 
not intersect F', it is possible to produce a subgraph that approximates the behavior of the original 
graph from the perspective of the given modulus problem. 
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